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Abstract Since their discovery, spicules have attracted increased attention as en- 
ergy/mass bridges between the dense and dynamic photosphere and the tenuous hot 
solar corona. Mechanical energy of photospheric random and coherent motions can be 
guided by magnetic field lines, spanning from the interior to the upper parts of the solar 
atmosphere, in the form of waves and oscillations. Since spicules are one of the most 
pronounced features of the chromosphere, the energy transport they participate in can 
be traced by the observations of their oscillatory motions. Oscillations in spicules have 
been observed for a long time. However the recent high-resolutions and high-cadence 
space and ground based facilities with superb spatial, temporal and spectral capacities 
brought new aspects in the research of spicule dynamics. Here we review the progress 
made in imaging and spectroscopic observations of waves and oscillations in spicules. 
The observations are accompanied by a discussion on theoretical modelling and inter- 
pretations of these oscillations. Finally, we embark on the recent developments made 
on the presence and role of Alfven and kink waves in spicules. We also address the ex- 
tensive debate made on the Alfven versus kink waves in the context of the explanation 
of the observed transverse oscillations of spicule axes. 
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1 Introduction 

The rapid rise of plasma temperature up to 1 MK from tlie solar photosphere towards 
the corona is still an unresolved problem in solar physics. It is clear that the mechanical 
energy of sub-photospheric motions is transported somehow into the corona, where it 
may be dissipated leading to the heating of the ambient plasma. A possible scenario of 
energy transport is that the convective motions and solar global oscillations may excite 
magnetohydrodynamic (MHD) waves in the photosphere, which may then propagate 
through the chromosphere carrying relevant energy into the corona. It is of great desire 
that the energy transport process(es) can be tracked by observational evidence of the 
oscillatory phenomena in the chromosphere. For a detailed discussion about MHD 
wave heating and he a ting d iagnostics in the solar atmosphere see the recent work by 
iTarovan and Erdelvil (|2009l ). 

Much of the radiation from the upper chromosphere originates in spicules, which are 
grass-like spiky features seen in chromospheric spectral lines at the solar limb ( see Fig 



1). Th ese abundant and spiky features in the chromos p here were discovered b vlSecchi 



18771) and were named "spicules" bv iRobertj (Il945h . lBecker3 (Il968l . Il972l ): ISterling 
(12000) dedicated excellent reviews to summarizing the observational and theoretical 
views about spicules at that time. Since these reviews, many observational reports of 
oscillatory phenomena in spicules appeared in the scientific literature. In particular, 
it is anticipated that signatures of the energy transport by MHD waves through the 
chromosphere may be detectable in the dynamics of spicules. A comprehensive review 
summarizing the current views about the observed waves and oscillations in spicules, to 
the best of our knowledge, is still lacking and such a summary has not been published 
yet in the literature. Here we aim to fill this gap. 

The goal of this review is to collect the reported observations about oscillations 
and waves in spicules, so that an interested reader could have a general view of the 
current standing of this problem. Here, we concentrate only on observed oscillatory and 
wave phenomena of spicules and their interpretations. We are not concerned about the 
models of spicule generation mechanisms; th e inte rested reader may find these la tter 
topics in the recent review bv lSterlind l|2000l ) or in lDe Pontieu and Erdelvil l|2006l ). 



Section 2 is a short summary about the general properties of spicules. Section 3 
describes the oscillation events reported so far for solar limb spicules. Section 4 outlines 
the views and discussions about the interpretation of spicule oscillations and Section 5 
summarizes the main results and suggests future directions of research. 



2 General properties of spicules 

Spicules appear as grass-like, thin and elongated structures in images of the solar lower 
atmosphere and they are usually detected in chromospheric Ha, D3 and Ca II H lines. 
These spiky dynamic jets are propelled upwards (at speeds of about 20 km s^^) from 
the solar surface (photosphere) into the magnetized low a tmosphere of the Sun. Accord- 
ing to early, but still valid estimates by .Withbroel ( 19831 ) spicules carry a mass flux of 



about two orders of magnitude that of the solar wind into the low solar corona. With 
diameters close to observati onal limits (<500 km ), sp icules have been l argely unex- 
plained. The suggestion by D e Pontieu et all ( 2004h and lDe Pontieu and Er delvi (2006j) 



of channeling photospheric motion, i.e. the superposition of solar global oscillations and 
convective turbulence, has opened new avenues in the interpretation of spicule dynamics 



Fig. 1 High resolution image of spicules at the solar limb in Ca II H line taken by Solar 
Optical Telescope (SOT) on board of Hinode spacecraft (November 22, 2006). 



(see a lso Hansteei^_et_^ 



( 20071 ) : [Heggland et al. 



] (l2006l ): lDe"Pontieu et aP (|2007ch : lRouppe van der Voort et all 
( 20071 )). The real strength of the observations and the forward 
modelling by De Pontieu and Erdelyi is that, as opposed to earlier existing models, 
they could account simultaneo usly for spicule ub i quity, evolution, energetics and the 
recently discovered periodicity (|De Pontieu et al ] ll2003al lbl)) of spicules. 

Excellent summaries about the general properties of spicule ( labelled t hese days as 
type I spicules) have been presented almost forty years ago bv lBeckerd (Il968l . 1 19721 ) 
and we broadly recall these findings here. Moreover, type II spicules were recently 
discovered with H i node a nd have very different properties from the classical spicules 
( De Pontieu et al.l l|2007bl )). 



2.1 Diameter 

Mea sured ra n ge of spicule diameter from ground based observations was ^ 700-2500 



Mea sured ra n ge ot spicule diameter trom ground based observations was ^ (\)0-Za 
km iBeckerd (|l968l )). The general view was that the diameter varies from spicule 



to 



spicule having the values from 400 km to 1500 km. The spicules seemed to be wider in 
Ca II H line than in Ha jBeckersI (Il972l )). However, the unprecedentedly high spatial 



resolution of Solar Optical Telescope on board of Hinode spacecraft (0.05 arc sec for 
Ca II H and 0.08 arc sec for Hq) revealed fine structure of spicules. Fig 1. shows 
this fine structure of spicules in Ca II H line at the solar limb taken by Hinode/SOT 
on November 11, 2006. The type II spicules disc overed by Hino de/SOT have smaller 



diameters (< 200 km) in Ca II H line (|Po Pont ieu et al.l 1 2007bl )). while the diameter 



of spicules in Ha line seems to be wider ~ 350-400 km (s ee also the high resoluti on 



image in Fig. [2] taken by Swedish Solar Telescope, courtesv lDe Pontieu et al] ( 200' 
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Fig. 2 High resolution image of spicules on solar disc taken by Swedish Solar Telescope (SST) 
in La Palma, adopted from .Dc Pontieu ct al. (.2004 ). 



2.2 Length 



The upper part of spicules continuously fade away with height, therefore the length is 
difficult to determine with precision. Generally, the top of a spicule is defined as the 
height where the spicule beco mes invi s ible. T he mean length of classical spicules varies 
from 5000 to 9000 km in Ha (jSeckersI j 19721 ')') and may reach to 7000-11000 km heights 
from the limb when observing by ground based coronagraphs. On the other hand, the 
type II spicules dominate in lower heights: they are tallest in coronal holes reaching 
heights of 5000 km or more, while in quiet Sun regions they reach lengths of ord er 
several megameters and they are shorter in active regions f'Pc Ponti eu et all (I2q07bh ). 
Additionally, very long spicules, called as macro spicules by Bochlin et al.l l|l975h with 



typical length of u p to 40 Mm are frequentl y observed mostly near t h e polar regions 
as reported by e.g. 'Pike and Harrison' ( 1997): IPike and Mason* 1*1998*): 'Paneriee et al.l 



2000); Parcnti ct al. (2002); Yaniauch i et al 
20061 ) : 1o^hea et all l|2005l ): IScullion et al.l 



2005 ) ; Doylo et al. (2 005); Madiarska et all 
20091 ): iNishizuka etlo] ijiooi) . 



2.3 Temperature and Densities 



Spicules have the temperatures and densities typical to the values of the chromospheric 
plasmas. Table [T] summarizes the typical electron temperatur es (Te) and num ber densi- 
ties (ne) of spicule values at different heights above the limb ( Beckers! ( 19681 )). Caution 
has to be exercised as values at 2000 and 10000 km heights are unreliable because of 
insufficient data. Typical electron densi t y at t he heights, where spicules are observed, is 
I AschwandenI 1 20041 ) . Fig 1.19 therein), therefore spicules are 



much lower 



10^ cm" 
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much denser than their surroundings. Matsuno and Hiravamal ( 19881 ) estimated lower 
temperatures (~ 5000-8000 K) than ISeckersI ( "gsd l: in general, spicule temperature 
seems to be much lower than that of the surrounding coronal plasma. 



Table 1 Electron temperatures and densities inferred from spicule emission, after iBeckerj 
l ll968f) 



h(km) 


Te{K) 


rie (cm ^) 


2 000 


17 000 


22x1010 


4 000 


17 000 


20xl0l0 


6 000 


14 000 


11.5x10"' 


8 000 


15 000 


6.5xl0i° 


10 000 


15 000 


3.5x1010 



2.4 Life time and motions 



The change of spicule length has been studied by many authors (see lBeckerd ( 19721 ') and 



references therein). The general opinion is that spicules rise upwards with an average 
speed of 20-25 km/s, reach the height 9000-10000 km, and then either fade or descend 
back to the photosphere with the same speed. The typical life time of classical spicules 
is 5-15 mins, but some spicules may live longer or shorter. On the other hand, the type 
II spicules from Hinode/S OT have much shorter lif e time, about 10-150 s and higher 
velocities of 50-100 km/s jPe Pontieu et all (l2007bl )'). Measurements of Doppler shifts 



in classical spicule spectra revealed the velocity of 25 km/s, similar to the apparent 
speed, therefore it was suggested that the apparent motion is real. However, it is 
also possible that observed Doppler shifts partly correspond to the transverse motions 
of the spicule axis and not to the actual movement along the axis. Such transversal 
motions can be caused due to the propagation of e.g. waves in spicules. These periodic 
perturbations are the subject of our discussion in the remaining part of the paper. 

It must be mentioned here, that some observations show the tilt of spicule spectra 
relative to the direction of dispersion, which was explained as due to the rotation of 



spicules around their axes IjPasachoff et al. I j 19681) : lRompoi^ (|l975l ): IPishkald (|l994l )) 



In this regards, it is interesting to note the recent SOHO/CDS observations, which 
also suggest the rotation in macro-spicules, interpreted as a sort of giant solar tornado 
( Pike and MasonI l|l998h ). 



3 Oscillations in solar limb spicules 

Oscillations in solar limb spicules can be detected either by imaging or spectroscopic ob- 
servations. Imaging observations may reveal the oscillations in spicule intensity and the 
visual periodic displacement of their axes. Imaging observations became especially im- 
portant after the recently launched Hinode spacecraft. SOT (Solar Optical Telescope) 
on board of Hinode gives unprecedented high spatial resolution images of chromosphere 
(see Fig. 1). On the other hand, the spectroscopic observations may give valuable in- 
formation about spicules through the variation of line profile. Variations in Doppler 
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shift of spectral lines can provide information about the line-of-sight velocity. Through 
spectral line broadening it is possible to estimate the non-thermal rotational veloci- 
ties leading to the i ndirec t observations of e.g. torsi onal Alfven wave s as suggested by 
lErdelvi and FedunI (|2007l ). and reported recently bv lJess et all (|2009l ) in the context of 
a flux tube connecting the photosphere and the chromosphere. Jess et al. used the tech- 
nique of analysing Doppler-shift variations of spectral lines, based on the optically thick 
Ha line for which a straightforward interpretation of linewidth changes and intensity 
changes in terms of velocity and density are sometimes very difficult and appropriate 
caution has to be exercised, and detected oscillatory phenomena associated with a large 
bright point group, located near solar disk centre. Wavelet analysis reveals full-width 
half-maximum oscillations with periodicities ranging from 126 to 700 s originating 
above the bright point, with significance levels exceeding 99%. These oscillations, 2.7 
km s~^ in amplitude, are coupled with chromospheric line-of-sight Doppler velocities 
with an average blue-shift of 23 km s^^. The lack of co-spatial intensity oscillations 
and transversal displacements rule out the presence of magneto-acoustic wave modes. 
The oscillations are interpreted as a signature of torsional Alfven waves, produced by 
a torsional twist of ±25 degrees. A phase shift of 180 degrees across the diameter of 
the bright point suggests these Alfven oscillations are induced globally throughout the 
entire brightening. The estimated energy flux associated with this wave mode seems 
to be sufficient for the heating of the solar corona, once dissipated. The question self- 
evidently emerges: Could spicules guide similar (torsional) Alfven waves and leak them 
to the upper solar atmosphere? 

Let us return to the possibility of intensity variations of spicules. Variation of line 
intensity indicates the propagation of compressible waves. And finally, the visible dis- 
placement of spicule axis may reveal the transverse waves and oscillations in spicules. 
Note that ground based coronagraphs can play an especially important role in spec- 
troscopic observations. Spatial resolutions of ground based coronagraphs reach to ~ 
1 arcsec, which is less than the resolution of Hinode/SOT. However, observations on 
coronagraphs usually are performed at 5000-10000 km above the surface, where spicules 
are less frequent and well separated. Hence, the appearance of two unresolved spicules 
inside 1 arcsec is unlikely (note, that the newly discovered type H spicules do not 
usually reach these heights). Therefore, the spectroscopic observations give valuable 
information about the dynamics of individual spicule at higher heights. 

In this section we briefiy outline almost all oscillation events in solar limb spicules 
reported so far in literature. Each observation is described in a separate subsection. 
At the end of the section, we summarize the information gathered about the typical 
observed periods in limb spicules. 



3.1 iNikolskv and Sazanovl (|l967l ') 



A very early, to the best of our knowledge the first, modern account of Dop pler shift 
temporal variation in spicules was reported bv INikolskv and Sazano 3 (| 19671 ). A set of 
spectrograms of chromospheric spicules in Ha line were obtained with the coronagraph 
of the Institute of Terrestrial Magnetism (Russia) on 1 August 1964. The Ha profiles 
and radial velocities of 11 different spicules were successfully derived from successive 
Ha line spectra formed at a height of ~6000 km above the solar limb. The time duration 
of observations was 8 mins with the intervals between exposures 10-40 s. Fig. [3]displays 
the time evolution of radial velocity and half width of Ha line for 10 different spicules. 
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Fig. 3 Temporal variations in radial veloc ities (solid curves) and half w idths (dotted curves) 
of 10 individual Ha spicules, adapted from [Nikolskv and Sazanovl ||1967D . 



Quasi-periodic oscillations are clearly seen. The authors concluded that the radial 
velocities vary randomly with time with a mean period of ~ 1 min. The amplitude 
of the oscillations are within 10-15 km/s. The half width of Hq line profile also tends 
to oscillate with a period similar to the mean period of ~ 1 min. The periodicity is 
similar to the time scale of type II spicules. Therefore, this may give the idea that the 
observed temporal variations are caused by the type II spicules. However, the type II 
spicules usually do not reach the observed heights. Therefore, the observed oscillations 
are unlikely to be connected to their activity. 



3.2 IPasachoff et all (|l968h 

IPasachoff et al.l ( 1968l l analyzed the high-dispersion spectra of the solar chromosphere 
obtained at the Sacramento Peak Observatory in several spectral regions separately 
during the summer of 1965. The observations were carried out simultaneously at two 
heights in the solar chromosphere separated by several thousands of kilometers. 

The time series of Doppler velocities in H-line of ionized Calcium, Ca II, for 10 
different spicules at the height of 5000 km are shown in Fig. [l] The exposure time for 
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Fig. 4 Temporal variation of Doppler velocities in H-line of Ca II from lPasachoff et al.l ||1968 



H-line was 13 s and the spatial resolution was less than 2 arc sec. lPasachoff et alj ( 19681 ) 
were searching for the sign reversal of Doppler velocities in order to determine the ris- 
ing/falling stages of spicule evolution. Indeed, some features show the sign reversal, but 
the common property is the clear quasi-periodic temporal variation of Doppler veloci- 
ties with periods of 3-7 min. The ampl itudes of oscillations ar e rather high, though still 
being within the range of 10-20 km/s. IPasachoff et alJ dlOGSl ) interpreted the detected 
temporal variation as motion along the spicule axis, but transverse oscillations also 
cannot be ruled out. 



3.3 IWeartI USM 



Observations have been carried out with the Mount Wilson Solar Tower Telescope 
during the period of 10 September - 13 October, 1967. Time sequences of Ha spectra 
with time lapse rates of 5 to 15 s have been obtained corresponding to height 5000-6000 
km height above the solar limb. 

The author reported that, both Doppler velocity and horizontal motion of spicules 
as a whole have significant input into spicule dynamics. In at least two cases, the author 
found that the combined motion indicate movement of a gas in an arc of a horizontal 
circle, firstly towards the observer, followed by sideways, finally away. Weart concluded 
that only true transverse motion could explain the observed pattern of motion. 

The power spectrum of temporal variations resembled the familiar 1/frequency 
curve, typical to many types of random motions. Substantial power was found to be 
concentrated at periods of 1, 2.5 and 10 minutes. However, no statistically significant 
peaks were observed. Therefore, it was concluded that spicules move horizontally at 
random. 



3.4 iNikolskv and Platoval l|l97lh 



Detailed spectroscopic observations were carried out on 3 April 1969 with the 53 cm 
Lyot coronagraph mounted at the High Altitude Astronomical Station near Kislovodsk 
(Russia). 38 Hq spectrograms of the chromosphere were obtained during about a 21 
minute observing campaign at the height of 4200 km. The time interval between suc- 
cessive frames varied from 14 to 100 s being on the average about 30 s. The spatial 
resolution of observations was ~ 1 arcsec. Two comparatively stable spicules, present 




in almost all frames, were chosen as reference ones. Then, the variation of other spicules 
along the limb with respect to these "bench-mark" spicules were determined. Fig.[5](left 
panel) shows the position of several spicules vs time with respect to the "bench-mark" 
spicules. There is evidence of oscillations of spicule position along the limb. 

The distribution of periods of spicule oscillations along the limb is shown on the 
right panel of Fig. [5] The most probable period lies between 50-70 s and the authors 
concluded that spicules undergo transversal oscillations with a period of ~ 1 min. The 
amplitude of oscillations was estimated to be about 10-15 km s""'^. The observations 
have been performed at the height, where the type II spicules may reach, therefore it 
is possible that observed temporal variations here are connected to their activity. 
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3.5 iKulidzanishvili and Nikolskvl (|l978l l 



T he observational material described i n the previous subsection was re-analyzed later 
bv IKulidzanishvili and Nikolskvl 1 19781 ) in order to search for a longer periodicity. The 
authors also searched for signatures of possible oscillations in Doppler velocity, line 
width and intensity, respectively. Fig. |6] shows the distribution of the spicule number 
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Fig. 6 Distribution of periods of Dopple r velocity V r, intensity W and full- width at half- 
maximum A\, adapted from IKulidzanishvili and Nikolskv, (.1978i) . 



vs the observed periods of oscillation in line-of-sight velocity, line width and intensity. 
About 70% of the observed periods of Doppler shift oscillations are within 3-7 min. 
The same per cent of observed periods lies within 4-9 min in the spicule intensity and 
80% of observed periods of line width oscillations are within 3-7 min. 



3.6 lGadzhiev and Nikolskvl l|l982l ) 



Observations were also carried out with the 53 cm Lyot coronagraph at Shemakha 
Astrophysical Observatory resulting in Ha time series corresponding to a height of 4 
Mm above the solar limb. 26 spectrograms have been taken over an 8-min interval 
which gives ~ 20 s between consecutive frames. A total number of 15 spicules were 
investigated in details. Gadzhiev and Nikolsky analysed variations in Doppler velocity 
as well as in the tangential velocity, i.e. reflecting the visible displacement of spicule 
axes along the solar limb. 

The authors found that the spicules oscillate with typical periods of 3-6 mins, both 
in line-of-sight and tangential directions. Fig. [7] shows the time variation of line-of- 
sight and tangential velocities in one of the spicules. The periodicity in both velocity 
components is clearly visible. Gadzhiev and Nikolsky also constructed the trajectories 
of spicule motion by putting together both velocity components. They concluded that 
spicules undergo a cyclic motion as a whole on an ellipse with an average period of 4 
mins. The average amplitude of this cyclic motion was 11 km s~^. 
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Fig. 7 Time variation of the radial and tan gential v elocities Vr, Vt and the modulus V of the 
velocity vector for a spicule, adapted from .Gadzhiev and Nikolskv (.1982) . 

3.7 iKulidzanishvili and Zhugzhdal (|l983l ) 

Another spectroscopic observations were carried out with the 53 cm Lyot coronagraph 
mounted at the Abastumani Astrophysical Observatory (Georgia). A 22 minutes long 
time sequence (the time cadence is not known), taken in the Ha hne, was analyzed 
for a total number of 25 spicules. The statistically significant period of oscillations in 
intensity, line width and line-of-sight velocity was found to be ~ 5 min. 



3.8 iHasan and Keill (|l984l ) 



Observations were also carried out with the Vacuum Tower Telescope of the National 
Solar Observatory (USA) in the autumn of 1982. A set of Ha spectra corresponding to 
five slit positions above the solar limb were recorded every 8 s in order to investigate 
the temporal variations of spicules. The spatial resolution of these observations was 
better than 2 arcsec. 

Hasan and Keil detected the temporal variations of the line-of-sight velocity at 
two different heights for two spicules. The fine time resolution allowed them to discern 
small amplitude fiuctuations with periods of about 2-3 mins. 



3.9 IPapushev and Salakhutdinovl (|l994l ) 



Observations were carried out with the 53 Lyot coronagraph of Sayan Observatory 
located near Irkutsk (Russia). The spectroscopic time series in different spectral lines 
varied from several minutes to hours with an excellent temporal resolution of 10-20 s. 
The spatial resolution of observations was better than 1 arcsec. The spectra were 
simultaneously registered at three different heights (5000 -8000 km above the limb) 
above the limb with a three-level image sheer. 
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Fig. 8 Temporal variations of Ha line profile parameters and the line-of-a ight velocity of 
spicules at different heights, adapted from lPapushev and Salakhutdinovl | |1994D . The solid line 
corresponds to the height of 5 Mm and the dashed line to 8 Mm above the solar limb, respec- 
tively. 



Temporal variations of Ha line profile parameters and line-of-sight velocity for one 
of spicules at two different heights are shown on FiglH] The quasi periodic fiuctuations 
are clearly seen. Papushev and Salakhutdinov found that the oscillation periods lay 
between 80-120 sec. 



3.10 IXia et all l|2005l ) 



Xia et al. analyzed the time series of EUV spicules in two polar coronal holes ob- 
tained by the SUMER (Solar Ultraviolet Measurements of Emitted Radiation) camera 
on-board the SOHO (SOlar and Heliospheric Observatory) spacecraft. The spatial res- 
olution of the observations was 1 arcsec and the exposure time for difi'erent data sets 
varied as 15, 30 and 60 s. Fig.|§]shows Dopplergrams and radiance map for the C III 977 
A line (left panel). The right panel shows the relative Doppler shifts at four different 
locations above the solar limb. The Doppler velocity and radiance indicate evidence of 
~5-min oscillations. 



3.11 iKukhianidze et all (|2006l ') 



Observations were carried out on 26 September 1981 with the 53 cm Lyot coronagraph 
of the Abastumani Astrophysical Observatory (the instrumental spectral resolution 
and dispersion in Ha are 0.04 A and 1 A/ mm correspondingly) at the solar limb. The 
scanning of height series began at the heigh t of 38 00 km measured from the photo- 
sphere, and continued upwards iKhutsishvilil l|l986l )). The chromospheric Ha line was 



used again to observe solar limb spicules at 8 different heights. The distance between 
neighbouring heights was l" (which was the spatial resolution of observations), thus 
the distance of ~3800-8700 km above the photosphere was covered. The exposure time 
was 0.4 s at four lower heights and 0.8 s at higher ones. The total time duration of each 
height series was 7 s. Con secutive height series began immediately, once a sequence was 
completed. Kukhianidze et al. (2006) analyzed the spatial distribution of Doppler ve- 
locities in selected Ha height series. Nearly 20% of the measured height series showed 
a periodic spatial distributions in the Doppler velocities. A typical Doppler velocity 
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Fig. 9 Time series Dopplergram and radiance map for t he C III 977 A line showing levels 
of radiance in logarithmic scale (adapted from IXia et al.l ||2005D ). The right panel shows the 
relative Doppler shift at four location above the limb. 
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Fig. 10 The Doppler velocity spatial distributions for one of the height series from 
iKukhianidze et al.l l|2006l ) . The marked dots indicate the observed heights. 



spatial distributions for one of the height series is shown in Fig. 1101 which shows a 
periodic behavior. The authors suggested that the spatial distribution was caused by 
transverse kink waves. The wavelength was estimated to be ~3.5 Mm. The period of 
waves was estimated to be in the range of 35-70 s. 
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3.12 IZaaarashvili et al.1 (|2007l l 



Zaqarashvili et al. analysed the same observational data obt ained by the 53 cm Lyot 
coronagraph of the Abastumani Astrophysical Observatory bv lKhutsish^diil l|l986l) . Us- 



ing the height series, they constructed continuous time series of Hq spectra with an 
interval of ~7-8 s between consecutive measurements at each height. The time series 
cover almost the entire lifespan (from 7 to 15 mins) of several spicules. Figure 1111 




Fig. 11 Left: Doppler vel ocity time series at th e heights of 5200 and 5900 km in one of 
the spicules, adapted from IZaqarashvili et al.l (12007). The time interval between consecutive 
measurements is ~8 s. Right: power spectra of Doppler velocity oscillations from the time 
series. The dotted lines in both plots show 95.5% and 98% confidence levels, respectively. 
There is the clear evidence of oscillations with periods of 180 and 30 s at both heights. 



shows the Doppler velocity time series at two different heights above the solar limb 
in one of the spicules (left panel) . The time series show the evidence of quasi-periodic 
oscillations in line-of-sight velocity. The power spectra resulted from Discrete Fourier 
Transform (DFT) analyses of the time series are presented in the right panel. The most 
pronounced periods at both heights are 180 and 30 s. The oscillation with the period 
of 90 s is also seen but preferably at higher heights (note the small peak at the lower 
height as well). 

The power spectra of Doppler velocity oscillations in two other spicules at the 
heights of 5200 km (lower panels) and 5900 km (upper panels) are plotted on Figure [T^ 
One of the spicules (left panels) shows the two clear oscillation periods of 120 and 80 
s at both heights. Both periods are above the 98% confidence level. Another spicule 
undergoes oscillations with periods ~110 and ~40 s, respectively. 

Zaqarashvili et al. also presented the results of DFT for 32 different time series as a 
histogram of all the oscillation periods above the 95.5% confidence level (see Figure [T3)). 
Almost half of the oscillatory periods are located in the period range of 18-55 s. Another 
interesting range of oscillatory periods is at 75-110 s, with a clear peak at the period 
of 90 s. Note that there is a further interesting period peak at 178 s as well, which is 
interpreted as a clear evidence of the well-known 3 min oscillations ubiquitous in the 
lower solar atmosphere. 
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Fig. 12 Power spectr a of Doppler veloci t y osci llations in another two spicules at the heights 
of 5200 and 5900 km l IZaQarastivili et al.l | |2007|) '). The dotted lines in both plots show 95.5% 
and 98% confidence levels. 




18 26 34 42 50 58 66 74 82 90 98 106 114 122 130 138 146 154 162 170 178 



periods (s) 

Fig. 13 Histogram of all osci llation periods that are above 95.5% confidence level, adapted 
from [Zagarashvili et al.l l l2007f) . The horizontal axis shows the oscillation periods in seconds, 
while the vertical axis shows the number of corresponding periods. 



3.13 iDe Pontieu et all (|2007al ) 



These authors analysed the time series of Ca II H-line images taken with SOT (Solar 
Optical Telescope) onboard the recently launched Hinode satellite. SOT has high spa- 
tial (0.2 arcsec) and temporal (5 s) resolutions, and is capable of providing extremely 
useful information about the dynamics of the lower atmosphere. 

lOe Pontieu et ahl (l2007al ) found that many of the chromospheric spicules (both, 
type I and II) undergo substantial transverse displacements, with amplitudes of the 
order of 500 to 1000 km during their relatively short lifetime of 10 to 300 s. They 
also reported that some longer-lived spicules undergo periodic motions in a direction 
perpendicular to their axes. Fig. 1141 shows an example of transverse d isplacement of a 
spicule with a period of ~3 mins. In general. lOe Pontieu et all (|2007al ) found it difflcult 
to determine the periods of transverse motions due to the short life time of spicules. 
Fig. 1151 shows the fine structure of spicules in Ca II H-line and the time-distance plot 
along the cut parallel to the solar surface. The plots reveal the predominant linear 
motion of spicules rather than oscillatory ones. However, comparing these excellent 
quality observations to their Monte-Carlo simulations of numerically modelling swaying 
spicule motion, they suggested that the most expected periods of transverse oscillations 
should lay between 100 and 500 s, which they interpreted as signature of Alfven waves. 
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Fig. 14 An example of the transversal motion of a spicule obtained with Hinode/SOT by 
iDe Pontieu et aU ll2007aD . Panel A shows the intensity as a function of time along the spatial 
cut indicated by a white line on panels B-G. Panels B-G represent the time sequence of Ca II 
H images. 




Fig. 15 Transverse motion of many spicules, from lDe Pontieu et al.l (|2007af ). Panel A shows 
an image of the Hinode/SOT Ca II chromospheric line. Panels B and D are time-distance plots 
along the cuts labeled by 1 and 2 on panel A. Panels C and E are similar cuts but reproduced 
by Monte-Carlo numerical simulations of spicule motions. 



3.14 Summary of observed oscillatory phenomena 

All the observations allow us to make a short summary of oscillatory phenomena de- 
tected in limb spicules. The results are gathered in Table [51 Oscillations are more 
frequently observed in the Doppler velocity, which points towards transversal motions 
as the observations are performed at the solar limb. The longitudinal velocity compo- 
nent may also take a part in shaping the Doppler velocity oscillations as spicules are 
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Table 2 Summary of observed oscillatory periods in solar limb spicules 







Dopplcr 


displacement 


intensity 


speed 


line 


Nikolskv and Sazanov ^^67) 




1 min 


1 min 






tia 


Pasachoff ct al. (1968) 




3-7 min 






> 90 km/s 


Ca II 


Wcart (1970) 




random 


random 






Hq 


Nikolskv and Platova (1971) 






50-70 s 






Ha 


j<Culidzanishvili and Nikols^k^ 1 


1978) 


3-7 min 




3-7 min 




Ha 


Gadzhicv and Nikolsky ( 1982' 




3-6 min 


3-6 min 






Ha 


j<ulidzanishvili and Zliugzhda 


(1983) 


5-min 




5-min 




Ha 


Hasan and Kcil (1984) 




2-3 min 






> 300 km/s 


Ha 


PaDushev and Salakhutdinov 


fl994) 


80-120 s 


80-120 


80-120 


> 300 km/s 


Ha 


Xia et ak (2005) 




5-min 




5-min 




EUV 


Kukhianidze et al. (2006) 




35-70 s 






~ 80 km/s 


Ha 


Zaaaraslivili et al. (2007) 




30-110 s 






~ 110 km/s 


Ha 


De Pontieu et al. (2007a) 






100-500 s 




50-200 km/s 


Ca II 



generally tilted away from the vertical. However, the transverse component seems to 
be the more important and determinant component in these oscillations as the visi- 
ble displacement along the limb is also frequently reported. The observed periods can 
be formally divided into two groups: those with shorter periods (<2-min) and those 
with longer periods > 2-min. The most frequently observed oscillations are within 
the period ranges of 3 — 7 min and 50 — 110 s. The two groups of oscillations are 
possibly caused by different physical mechanisms, but this issue needs further studies 
before one can conclude. Intensity oscillations are observed mostly with ~5-min pe- 
riod, which may indicate t h eir connection to the global photosphe ric 5-min oscillations 
( De Pontieu et al.l (|2003al lbl. [20o3 ): lDe Pontieu and Erdelvil (|2006h ). However, ~ 5-min 
oscillations in Doppler velocity and visible displacement can be caused by transverse 
waves (kink or Alfven). Oscillations with other periods seem to be connected mostly 
with transverse displacement of spicule axes, which can be caused by kink or Alfven 
waves. These possibilities are discussed later in details. 

In order to infer more precise information about the oscillatory motions in the lower 
solar atmosphere, in particular in chromospheric spicules, it is of vital importance to 
determine whether the oscillations are caused by propagating or standing wave pat- 
terns. Unfortunately, this is a very challenging task from an observational perspective, 
and it is a difficult task from a theorist's point of view as well. To the best of our 
knowledge, very little is known about th ese mechanism s appl icable to the highly strat- 
ified spicule geometry in the literature. lErdelvi et al. I (|2007l ) give an insight into the 
theoretical and numerical difficulties of studying waves and oscillations in the highly 
stratified lower solar atmosphere. In spite of these obstacles, some conclusions still can 
be drawn and we discuss them in the next section. 



4 Propagation 

The propagation of disturbances along spicules can be deduced when observations 
are performed at least at two different heights. Then the propagation speed can be 
estimated from the phase difference between oscillations at different heights (see Table 
II). 
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Several authors rep orted the possible propagation speeds of disturbances in spicules. 
IPasachoff et alJ ( 1968l l found that velocity changes occur simultaneously, to within 20 s, 
at two distinct heights separated by 1800 km. They concluded that the propagation 
velocities should be more than 90 km s"^ iHasan and Keill (|l984l ') suggested the prop- 
agation of signals from lower to higher heights with a n estimated velocity of more 
than 300 km s~^. IPapushev and Salakhutdinol (| 19941 ) studied the phase delays of 
fluctuations at different heights and also concluded that the propagation speeds should 
exceed 300 km s~^. However, it must be noted that a standing oscillation pattern may 
also be responsible to ex plain these phenomena. Indeed, Hinode/SOT observations 
( De Pontieu et al.l (|2007al 'l) show some evidence of upward and downward propaga ting 
waves, with some partially standing waves being observed. De P ontieu ct al. (2 007al ) 
estimated wave phase speed in the range of 50-200 km/s. 

Detailed reports abo ut wave prop agation were presented by e.g. iKukhianidze et al.l 
(I2OO6) and Zaaarashvili et al.l (2003) through analyzing the consecutive height series. 
^ukhianidze et al.1 i20od)^e;;^;rthree consecuti;e height series of Doppler velocities 
in a spicule, which show that the maximum of the Doppler velocity moves up in time 
(see Fig. I16p . The authors suggested that this may indicate a wave phase propagation. 




10 20 10 20 10 15 20 
velocity (km/s) velocity (km/s) velocity (km/s) 



Fig. 16 Three consecut i ve hei ght series of Doppler velocities in one of the spicules, adapted 
from IKukhianidze et all (120061) . The time difference between the consecutive plots is ^8 s. 
The maximum of Doppler velocity moves up in consecutive height series, which most probably 
indicates wave propagation. 



The phase is displaced at ~1500 km in about 18 s giving the phase speed of ~80 km s""'^, 
ver y comparable to the expe cted kink or Alfven speed at these heights. 

IZaqarashvili et all (2003) presented a Fourier power as a function of frequency and 
heights for two different spicules shown on Fig. [17] (left panel). There is clear evidence of 
persisting oscillations along the full length of both spicules. The plot of the first spicule 
shows the long white feature (feature A) located just above the frequency 0.01 s~^. This 
is the oscillation with the period of ~80 s and it persists along the spicule. The most 
pronounced feature (feature B) in the plot of the second spicule is colour-code indicated 
by a long brighter trend located just above the frequency of 0.02 s~^ and persisting 
along almost the whole spicule. This is the oscillation with a period of 44 s. Then, 
Zaqarashvili et al. calculated the relative Fourier phase between heights for the most 
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leigfj Heights 

Fig. 17 Left; Fourier power expressed i n confidence levels a s function of frequency and heights 
for two different spicules, adapted fromEag arashvili et alj (2007). Brighter points correspond 
to higher power, and darker points correspond to lower one. The label 1 on the power scale 
(right plots) corresponds to the level of 100% confidence. Right: Relative Fourier phase as a 
function of height for oscillations (a) with ^80 s period in the first spicule (feature A) and (b) 
44 s period in the second spicule (feature B). The distance between consecutive heights is ~ 1 
arc sec. 



pronounced features (features A and B). Right panel of Figure [17] shows the relative 
Fourier phase as a function of heights for (a) feature A and (b) feature B. In spite of 
the apparent linear behaviour of the phase difference, there is practically almost no 
phase difference between oscillations at different heights for feature A, which probably 
indicates a standing-wave like pattern with period of ~ 80-90 s. On the contrary, there 
is a significant linear phase shift in plot (b), which indicates a propagating pattern 
with an estimated period of ^40-45 s. Therefore, the authors concluded that the first 
spicule shows the standing-wave like pattern (or a wave propagation almost along line 
of sight, which seems unlikely though cannot be ruled out) while the second spicule 
shows a propagating wave pattern. Zaqarashvili et al. estimated the propagation speed 
as ~110 km s~^. 

Hence, based on these few observations we conclude that the propagation speed 
of disturbances in solar limb spicules is quiet high and may exceed 100 km . This 
indicates that the disturbances are of magnetic origin and they propagate with chro- 
mospheric Alfven or kink speeds that exceeds the local sound speed in spicules at these 
heights (but, note that the sound speed outside spicules has almost coronal values). 



5 Possible interpretation of observed oscillations 

It is clear that the observed transverse oscillations of spicule axis can be explained and 
interpreted by the waves propagating along the spicule. The two type of waves respon- 
sible for periodic transverse displacement of spicule axis are: kink or Alfvon waves. 
If spicules are modelled as plasma jets being shot along magnetic fiux tube, the n th e 
transverse oscillati o ns co ul d be caused by MHD kink waves Ikukhianidze et all hooj ): 
IZagarashvili et ah I l|2007l ): lErdelvi and FedunI (|2007l ): lAiabshirizadeh et all (|2009l )). If 
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a spicule is not stable w ave guide for the tube w aves, then the oscillations can be 
caused by Alfven waves ( De Pontieu et aP 1 2007ah ). Before we embark on the inter- 
pretation of spicule oscillations, let us briefly summarise the possible MHD modes 
and their main properties in an inhomogeneous magnetised plasma under lower solar 
atmospheric conditions. 



5.1 MHD waves in a uniform magnetic cylinder in the lower solar atmosphere 




Fig. 18 Left: Magnetic flux tube showing a snapshot of Alfven wave perturbation propagat- 
ing in the longitudinal z-direction along fleld lines at the tube boundary. At a given height 
the Alfenic perturbations are torsional oscillations, i.e. oscillations are in the ip-direction, per- 
pendicular to the background field. Note that on the other hand MHD kink waves would force 
the tube axis to oscillate. Right: Snapshot showing Alfven waves propagating along a mag- 
netic discontinuity. Again, the key feature to note is that Alfvenic perturbations are within 
the magnetic surface (j/z-plane) at the discontinuity, perpendicular to the background field {y- 
direction), while the waves themselves propagate along the field lines (z-di rection). The MHD 
kink waves oscillate in xz-plane in this geometry. Image adapted from lErdelvi and FedunI 
l|2007D . 



To simplify the bewildering complexity of the dynamic solar a tmosphere, the con- 
cept o f magnetic flux tubes is often used. In a pioneering work by lid win and Roberta 
( 19831 ). using cylindrical coordinates, it was derived the dispersion relations of MHD 
waves propagating in cylindrical magnetic flux tubes. The main obstacle to be over- 
come when introducing the concept of flux tubes is the conversion from Cartesian to 
cylindrical coordinates. This change results involving Bessel functions in the dispersion 
relation which are not yet possible to be solved analytically without simplification, e.g. 
through incompressibility or long and short wav elength approximations. Let us sum- 
marise here the key steps of lEdwin and RobertsI (,^1983). Consider a uniform magnetic 
cylinder of magnetic field Bqz confined to a region of radius a, surrounded by a uni- 
form magnetic field SeZ (see Figure [TSk ) . To simplify the MHD equations we assume 
zero gravity, there are no dissipative effects and all the disturbances are linear and 
isentropic. Pressure (plasma and magnetic) balance at the boundary implies that 



Bq bI 

2^0 2^0 



(1) 
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where po and pe are the pressures inside and outside the tube. 

Linear perturbations about this equihbrium give the following pair of equations valid 
inside the tube, 

-i£)^-0, (3) 

where is the Laplacian operator in cylindrical coordinates (r, 0, z) and 

A = divv, r — z ■ curlv (4) 

for velocity v = {vr,v^, Vz). A similar pair of equations to ((2)1 and (|3]) are valid outside 
the tube. Fourier analysing we let 

A = R{r)exp[i(Lot + n(t> + kz). (5) 
Then equations ([2} and ([3]) give Bessel's equation satisfied by R{r) as follows 

d^R IdR . . 
where 

We have used the notation Vj^ for the Alfven speed, cp for the sound speed and ct 
for the characteristic tube speed (sub-Alfvenic), where ct = cova/{co + '^a)"^^^- To 
obtain a solution to (O bounded at the axis (r — 0) we must take 

[ Jn{nor), Uq = -rriQ > J 

where is an arbitrary constant and In, Jn are Bessel functions, see e.g. lAbramowitz and StegunI 
(1963), of order n. For a mode locked to the waveguide it is required that no energy 




propagates to or from the cylinder in the external region, i.e. the waves are evanescent 
outside the flux tube. Therefore we take 

R{r) = AiKnimer), r > a, (9) 

where Ai is a constant and 



which is taken to be positive (no leaky waves). Here, Ce stands for the sound speed 
outside the tube. Since we must have continuity of velocity component Vr and total 
pressure at the cylinder boundary r = a, this yields the dispersion relations 



/; 2 2 2\/i 2 2 2\ 
[k Ce - to ){k Vj^e ~ ^ 
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pQ(k - uj )me 



' KnijneO) 



/, 2 2 2n In^m-QO) 



for surface waves (ttiq > 0) and 



., 2 2 2n K'niniea) 2 2 2n ^Un-oa) 

Kn{mea) Jn{nQa) 



(11) 



(12) 



for body waves {niQ = — jiq < 0). The axisymmetric sausage mode is given by n = 0, 
while the well-observed kink mode (non-axisymmetric) is given by n = 1. Modes with 
n > 1 are called flute modes. Although the dispersion relations (|11[) and (|12p are 
complicated, finding the phase speed for e.g. kink waves with photospheric parameters 
simplifies matters considerably. 

Under lower solar atmospheric conditions, characterised by Ce > > co, possibly 
representative for spicules, sunspots or pores both the slow and fast bands have surface 
and body modes, respectively. The slow MHD waves are in a narrow band since cq « cy. 
The slow body waves are almost non-dispersive, whereas the almost identical slow 
surface sausage and kink modes are weakly dispersive (bottom zoomed out panel in 
Fig. I19p . On the other hand, if the characteristic speeds of a lower solar atmospheric 
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Fig. 19 The solution of the dispersion relations I I11I I- II121 in terms of phase speed {lo/k) of 
modes under photospheric conditions Ce > va > cq (all sp eeds are in km/s). The slow band is 
zoomed (lower panel). Image adapted from lErdelvil (|2008f l. 



fiux tube render as vj^ > Ce > cq, the fast body modes do not exist since they would 
become leaky mode solutions of the dispersion relations (|11|) - H12|) . The MHD modes 
propagating along the fiux tube in this case are shown in Fig. 1201 Note, there is little 
evidence, whether fiux tubes modelling spicules, are characterised by Ce > va > cq or 

VA > Ce > Cq. 
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Fig. 20 Similar to Fig. [T9]but for va> c^> cq- Image adapted from lErdelvil ll2008f) . 



In the incompressible limit, suitable for the description of kink waves and oscil- 
lations in their linear limit, (cg — > oo,Ce cxo),mo and me become |fc|. The kink 
and sausage modes are then given explicitly after some algebra. It is noted that the 
phase speed for the kink mode is not monotonic as a function of k but has a max- 
imum/minimum and the sausage mode is monotonically increasing/decreasing. This 
max/min feature of the kink mode is absent in the slab case, so it can be deduced 
to be a reflection of the geometry of the magnetic field. On the other hand, the in- 
compressible Alfven waves in magnetic tubes are polarized in the direction and do 
not lead to the displacement of the tube axes (see Figure ITSb). The propaga t ion o f 
tors ional Alfven waves in vertical magnetic tubes have been studied bv lHollwea ( 198ll ) 
and iHollweg et all (|l982| ) . 

Particularly interesting processes may take place near the heights, where the Alfven 

and sound waves have similar values of phase speed. This area with /? = !(/?= 

87rpn/j3n)i m ay locate somewhere between the photosphere and chromosphere. ISogdan et al.l 
(|2002| . I2OO3I ') have performed 2D numerical simulations in isothermal atmosphere and 



have s hown the coupling between different MHD wave modes near this area. Zagar ashvili and Robert^ 
( 20061 ) also show that the nonlinear coupling between Alfven and sound waves may take 
place there. 

Last but not least, for completeness and for the benefit of interested readers, we 
note that kink MHD wave propagation under solar coronal conditions is discussed in 
details bv lRuderman and Erdelvir(|2009l ) in this Volume. 



5.2 MHD kink waves 



Equipped with a clear understanding of the differences between kink and Alfvenic 
perturbations, as described above, let us now return to plausible interpretations of pe- 
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riodic spicular motions. Spicules are much denser than the surrounding medium above 
the height of 2000 km (see section 2.2), therefore they may be considered as magnetic 
tubes and the MHD wave theory can be applied in some format as outlined in the 
Sec. 15.11 Then, the observed periodic transverse displacement of the axis probably is 
due to the propagation of kink waves (see Figs. I18l andl2i p . Transverse kink w aves can 
be gen e rated in phot o spheric magnetic t ubes by buffeting of granular motions ( RobertsI 



be gen e rated m phot o sphenc magnetic t ubes by buitetmg ot granular motions l|Kobertsl 
( 19791 ') ; ISpruij (|l98l[ ) ; lOsin et al.l (1999') ) . The waves may then prop agate thro ugh th e 



stratified chromosphere (see, e.g. Hargreave^ ll2008l): iHargreaves and Erdelvii (120091)') 
and le ad to the observed oscillations ( Kukhianidze et al.l ( 2006 h_ Zagarashvili et al.l 
(2003)). The recent theoretical paper bv lAiabshirizadeh et al.l j2009[ ) also shows that 



the kink waves with periods of ~ 80-120 s may propagate in spicules at higher heights. 
Indeed, t he oscillations of Doppler velocity in spicules with similar periods were re- 
ported bv lZaaarashvili et al.l (20o3). If the velocity of the kink wave is polarized in the 



jNikolskv and Sazanov ( 1967l'):|Pasachoff et al.Nl968l'): Kuhdzanishvih and Nikolskvl (119781' 


ICadzhiev and Nikolskv (1982): Kuhdzanishvih and Zhugzhda 


(Il983t): Hasan and Keill 


(1984): PaDushev and Salakhutdinovl( 1994h: Xia et all (20051): 


Kukhianidze et all (20061): 



However, if the velocity is polarized in the per pendicular 



iZagarashvili et al.l 

plane t hen it results the visible displ a cement of spicule axis along the limb (iNikolskv and Sazanov 
19671 ) : lNikolskv and Platoval (|l97lh : lGadzhiev and Nikolskvl (|l982l ): |Papush ev and Salakhutdinovl 
19941 )). Kink wave propagation along a vertical thin magnetic flux tube embedded in 



observer 




spicule 



chromosphere 



Fig. 21 Schemat ic gicture of propagating kink waves in spicules, adapted from 

iKukhianidze et al.l ||2006|) . 



the stratified field-frcc atmosphere is governed by the Kloin-Gordon equation ('Rae and Robert: 
(1982); Spru it and Roberts (jl983); Hasan and Kalkofcn (1999); Zaqarashvili and Skhir tladzei 
( 20081 ) : lHargreavesi (|200j ): IHargreaves and Erdelvii (|2009l ): the latter authors have even 
considered a dissipative medium resulting in a governing equation of the type of Klein- 
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Gordon-Burgers equation) 



k k 



where Q = ^{z, t) exp (— z/4yl), cj, = Bq/ 47r(po + Pe) is the kink speed, yl is the den- 
sity scale height and = cj,/4yl is the gravitational cut-off frequency for isothermal 
atmosphere (temperature inside and outside the tube is assumed to be the same and 
homogeneous). Here £,{z,t) is the transversal displacement of the tube, Bq[z) is the 
tube magnetic field, pa{z) and pe{z) axe the plasma densities inside and outside the 
tube respectively (the magnetic field and densities are functions of z, while the kink 
speed Cfc is constant in the isothermal atmosphere). 

Eq. (I13|l yields simple harmonic solutions exp[i(a;t ± kzz)] with the dispersion re- 
lation 



^'^ - ^1 = clkl, (14) 



where uj is the wave frequency and kz is the wave number. The dispersion relation 
shows that waves with higher frequency than i7fc may propagate in the tube, while the 
lower frequency waves are evanescent. 




Kink waves cause the transverse displacement of the entire tube. The displace- 
ment of tube in a simple harmonic kink wave is polarized arbitrarily and the po- 
larization plane depends on the excitation source. The superposition of two or more 
kink waves polarized in different planes may give rise to the complex motion of the 
tube. The process is similar to the superposition of two plane electromagnetic waves, 
where the waves with the same amplitudes lead to the circular polarization, while 
the waves with different amplitudes lead to the elliptical polarization. Consider, for 
example, two harmonic kink waves with the same frequency but polarized in the xz 
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and yz planes: Ax = AxQCOsiut + kzz) and Ay — AyQsm{Lot + kzz). The superpo- 
sition of these waves sets up helical waves with a circular polarization if A^o = A^q 
( Zagarashvili and Skhirtladzel ( 
vertical, while the displacement remains constant (Fig. [22]). If Axo 7^ AyQ then the re- 
sulting wave is elliptically polarized. The superposition of few harmonics with different 
frequencies and polarizations may lead to an even more complex motion of the tube 
axis. 



20081 )). As a result, the tube axis rotates around the 



5.3 Alfven waves 



IPe Pontieu et ah! ( 2007a') suggested that the transverse displacement of spicule axis can 
be explained by the propagation of Alfven waves excited in the photosphere by granular 
motions or global oscillation patterns. They performed self-consistent 3D radiative 
MHD simulations ranging in the vertical direction from the convection zone up to 
the corona. A snapshot from the simulations is presented on the Fig. [23] (panel A). 
Their analysis shows that the field lines (red lines) in the corona, transition region, 
and chromosphere are continuously shaken and carry Alfven waves. Panels B and 
C are time-distance plots from the simulations and observations, respectively. From 
the comparison between simulations and observations the authors concluded that the 
period of Alfven waves should lay between 100 and 500 s. However, they suggested that 
very long-lived macro spicules show some eviden ce of Alfven waves with longer periods 
between 300 and 600 s. jPe Pontieu et al.1 l|2007ah also claim that their observations and 



simulations do not show the spicules as stable wave guide for kink waves. Therefore they 
argue that volume-filling Alfven waves cause the swaying of magnetic field lines back 
an d fo rth leading to the visible displacement of spicule axis. These claims are debated 
by lErd clvi and Fcdun (2007): "However, these observations also raise concerns about 
the applicability of the classical concept of a magnetic flux tube in the apparently 
very dynamic solar atmosphere, where these sliding jets were captured. In a classical 
magnetic flux tube, propagating Alfven waves along the tube would cause torsional 
oscillations (see Fig. I18h in this paper earlier). In this scenario, the only observational 
signature of Alfven waves would be spectral line broadening. Hinode/SOT does not 
have the appropriate instrumentation to carry out line width measurement^ On the 
other hand, if these cla ssical flux tubes did indeed exist, then the observations of 
iDe Pontieu et al. I l|2007al ) would be interpreted as kink waves (i.e., waves that displace 



the axis of symmetry of the flux tube like an S-shape) . More detailed observations are 
needed, perhaps jointly with STEREcQ, so that a full three-dimensional picture of 
wave propagation would emerge." 



5.4 Kink vs Alfven waves 



It is important to determine accurately which type of waves are responsible for t he 
transverse displacement of the tube axis (|Erdelv i and Fedunl (|2007l '): |jess et al.1 (|2009l )). 



^ We need to clarify here that, of course, Hinode/SOT can measure line width, however, not 
with the desired resolution and in the desired wavelength in the context of spicule oscillations, 
as opposed to e.g. the Rapid Oscillations in the Solar Atmosphere (ROSA) instrument mounted 
on the Dunn Solar Telescope, NSO, Sacramento Peak. 

^ Although STEREO has a limited spatial capacity for spicule observations, it can give a 
global view of the related MHD waveguide. 



27 





The existence and main features of wave modes depend on the properties of the medium 
where the waves propagate in (see Sec. I5.1|l . If spicules represents a magnetic flux 
tube, then MHD wave theory may permit the propagation of kink and torsional Alfven 
waves. The kink waves are global tube waves and cause the displacement of the tube 
as a whole. However, most importantly let us emphasise once more that torsiona l 
Alfven waves do not lead to th e displacement of t ube axis (jErdelvi and FedunI (|2007l ): 
IVan Do orssela ere et al.l ( 2008 al): |jess et all (|2009l )). Therefore, torsional Alfven waves 
are unlikely to be the reason for the observed oscillations. However, global Alfven waves, 
which may indeed fill a significant volume in and around spicules, may cause the global 
transverse oscillations of magnetic field lines. In this latter scenario spicules may Just 
simply follow the oscillations of field lines and this is what lOe Pontieu et ah 1 l|2007ah 
suggest. However, it must be mentioned here, that spicules may respond to the oscil- 
lation of magnetic field lines slower than the surrounding plasma d ue to their higher 
density i.e. larger inertia. This point is not taken into account by lOe Pontieu et aU 
(j20074 ) and it needs to be addressed in future models. 

However, there are two possible difficulties associated with the Alfven wave scenario. 
Firstly, if the oscillations of spicules are due to global Alfven perturbations, then the 
neighbouring spicules should show a coherent oscillation. However, Hinode movies in- 
dicate the opposite, i.e. spicule movements are random and there is no sign of coherent 
oscillations. This incoherence partly can be caused due to the fact that spicules seen 
at the limb are located in different parts along the line of sight. However, detailed 
analysis still can uncover the coherent oscillations and this point needs to be addressed 
in future. Secondly, it is not yet clear how the volume-filling Alfven waves would be 
generated in the photosphere, where the magnetic field is rooted and concentrated in 
fiux tubes. In this regards, it seems to be more plausible that the transverse pertur- 
bations propagate upwards in the form of kink waves, which may be transformed into 
Alfven-like waves in the chromosphere, where the magnetic field rapidly expands. The 
Alfven waves a lso can b e generated near P — 1 regi on due to various wave coupling 
processes (Bo gdan et all (^2002, 2003); Zaqarashvih a nd Roberts! (|2006l ')'). 
The kink wave scenario also has its own difficulties as the photospheric fiux tubes 
may be expanded in the chromosphere as we already noted above. This expansion may 
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cause certain difficulties for the kink wave to propagate. However, plasma /3 becomes 
less than unity in higher heights, which means that the concept of magnetic tube is 
changing compared to the photospheric conditions. Now, the magnetic tube means the 
higher concentration of density. This is exactly the case for spicules, as their density can 
be up to two order magnitude higher, than in the surrounding coronal plasma. There- 
fore, the spicules can be wave gui de of the kink waves (at least, in classical spicules) 
although lOe Pontieu et al] 1 2007ah claim opposite. Propagatio n of transverse pulse in 
the wa ve guide of enhanced density was recently studied by IVan Doorsselaere et al.l 
( 2008bl ). They show that the slab with enhanced density essentially trap the initial 
transverse pulse. It must be mentioned that the rapid disappearance of type II spicules 
in Ca II of Hinode/SOT does not immediately mean that spicules are not stable wave 
guides. The disappearance can be caused due to the increase of temperature inside 
spicule. But the tube (i.e. higher concentration of density) still may remain. In this 
regards, lErdelvi and FedunI (|2007l ). in the context of prominence oscillations, showed 
that kink waves can be easily guided. For the recently developed theory of transversal 
waves and oscillations in gravitationally and magnetically stratified flux tubes see, e.g. 
IVerth and Erdelvil l|2008h"lRuderman et al.l l|200d ): lAndries et al.l (|2009l ). The problem 
whether the spicule displacements are Alfvenic or kink waves is currently under de- 
bate and more sophisticated observat ions complemented by nu merical investigations 
are needed for a satisfactory solution l Erdelvi and FedunI 1 200?! )). 



5.5 Transverse pulse 



It must be recognised that simple harmonic waves can hardly be excited in the dynamic 
solar photosphere. A more realistic process of wave excitation is the impulsive bufi'eting 
of granules on an anchored magnetic flux tube, which may easily generate transverse 
pulses. Such pulses may propagate upwards in the strati fied atmosphere and leave the 
"wake" oscillating at cut -off frequency of kink waves ( Zagarashvili and Skhirtladzel 
(|2008l )). The wake may be also responsible for the observed transverse oscillations of 
spicule axes. 

For the sake of simplicity, let us consider _the simplest impulsive forcing in both 
time and coordinate. Then Eq. (|13|) looks as ffZagarashvili and Skhirtladza ()2008il) ) 



1 d^Q 



9t2 



-fQ = -AaS{mz), 



(15) 



where z > —oo, t > 0, Aq is a constant and the pulse is set off at t = 0, z = 0. 

The solution of this equation can be written, after e.g. (|Morse and FeshbachI ( 195^ )). 
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(16) 



where Jq and H are Bessel and Heaviside functions, respectively. Eq. (|16|l shows that 
the wave front propagates with the kink speed cj. (the first term), while the wake 
oscillating at the cut-off frequency Qj- is formed behind the wave front (the sec ond 
term ) and it decays as time progresses jRae and Roberts 1982; Sp ruit and Ro bertsI 
Ii983| :Ih asan and Kalkofenll 19991 : lHargreave3 l boOSljHargreaves and Erdelvl2009l ). Note, 
the actual mathematical form of the governing equation and its solution for kink 
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Fig. 24 Propagation of transverse kink p ulse in stratified vertical magnetic flux tubes, adapted 
from IZagarashvili and Skhirtladzd l l2008t) . A rapid propagation of the pulse is seen, which is 
followed by the oscillating wake. The pulse propagates with the kink speed and the wake 
oscillates with the cut-off period. The amplitude of the pulse (and wake) increases with height 
due to the decreasing density. 



and longitudinal oscillation s in a gravitationally st ratified and anchored magnetic flu x 
tube is very s imilar , see e.g. lSutmann et al. I (| 19981 ) : lMusielak and Ulmschneiderl (|200ll ): 
iBallai et al.1 l|2006l ) for more details. 

Fig. 1241 shows the plot of the transverse displacement (,{z,t) — Q{z,t) exp{z/4:A), 
where Q is expressed by the second term of Eq. (I16|) . A rapid propagation of the 
pulse is found, which is followed by the oscillating wake (the time is normalized 
by the cut-off period Tj, = 2n/f2i^). Just after the propagation of the pulse, the 
tube begins to oscillate with the cut-off period at each height. The amplitudes of 
pulse and wake increase upwards due to the density reduction, but the oscillations 
at each height decay i n tim e. A very similar an d resembling behaviour wa s foun d 
by iFleck an d Deubner' (1989'); ' Fleck and Schmitj l|l99lh ; ISchmitz and Fleckl (Il998l '): 
lErdelvi et al.. (,2007 ): Malins and Erdelvl (2007,) for longitudinal oscillations. 

Hence, the transverse and impulsive action on the magnetic tube at i = near 
the base of the photosphere (as set at 2 = 0) excites the upward propagating kink 
pulse, while the tube in the photosphere oscillates at the photospheric kink cut-off 
period, which can be estimated as ~8 min, using the photospheric scale height of 125 
km. Hence, the magnetic tube will oscillate with ~8 min period in the photosphere. 
When the pulse penetrates into the higher chromosphere, where the cut-off period is 
changed, this will influence the propagation of the pulse. Spicules with higher density 
concentrations than the ambient plasma may guide kink waves with the phase speed of 
25 km/s, which yields the cut-off period of about 500 s. Note, the Alfven cut-off period 
can be as short as 250 s. Therefore, the transverse pulse may set up the oscillating 
wake in the chromosphere with the period of 250-500 s. 
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It is most likely, that the anchored magnetic flux tube undergoes granular bufltet- 
ing from many directions. Therefore, the superposition of consecutive pu lses may set 
up a h elical motion of tube axis with the cut-off period ( Zaqarashvili and Skhirt ladzj 
2008 )). The helical motion of spicule axis first has been observed by Gadzhiev and Nikolskvl 



19821 ) from simultaneous observation of Doppler velocity and visual displacements. The 



recent Hinode/SOT movies also show complex motions of spicule axis. Fig. 1251 (upper 
panels) show the observed trajectories of spicule axis with respect to the photosphere, 
adopted from lCadzhiev and Nikolskvl 1 19821 '). The lower panel shows the superposition 
of two wakes at the height of 250 km above the photosphere, numerically simulated by 
IZaqarashvili and Skhirtladzel ( 2008h . by solving the Klein-Gordon governing equation. 
The first wake corresponds to the pulse imposed along the a;-direction and the second 
wake is the result of another pulse generated in the j/-direction with a different ampli- 
tude. The observed period of helical motions is in the range of 3-6 mins, which may 
correspond well with the kink cut-off period in the chro mosphere. A very similar phe- 
nomenon was recently observed in fibril orientation bv iKoza et al.l ( 20071 ) . Therefore, 





Fig. 25 Up per panels: Observed trajectories of the motion of two distinct spicules, adapted 
from lGadzhi ev and Nikolskv (1982). Lower panel: Simulated helical motion of the tube axis at 
the height of 250 km measured from the solar surface, due to the p ropagation of two consecutive 
transv erse pulses polarized in perpendicular planes, adapted from lZagarashvili and Skhirtladzg 

Hooi). 



the oscillation of wakes behind a tran sverse pulse may explain t he visible transverse dis- 
placement of spicule axis observed bv lDe Pontieu et all (l2007al ') with Hinode/SOT. The 
cut-off period is similar or slightly shorter than the mean life time of spicules, therefore 
the oscillations are difficult to detect as it is noted bv IPe Pontieu et^l . ( 2007a). It must 
be mentioned, however, that the photosphere/chromosphere is in much more complex 
and dynamic state than it is described by this simple theoretical approach. Therefore, 
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it is desirable to perform more sophisticated numerical simulations of transverse pulse 
propagation from the photosphere up to the corona. 

6 Summary of main results and targets for future work 

Spicules are one of the most plausible tracers and trackers of the energy coupling and 
energy transport from the lower solar photosphere towards the upper corona by means 
of MHD waves. These waves may induce the oscillatory phenomena in the chromo- 
sphere, which are frequently detected in limb spicules. Periodic perturbations, e.g. in 
forms of oscillations, are observed by both spectroscopic and imaging observations. Let 
us summarize the main observed oscillatory phenomena (see also Table [2}: 

— Oscillations in limb spicules are more frequently observed in Doppler shifts and in 
the visible displacement of spicule axis, which probably indicate the presence of 
transversal motions of spicules as a whole. 

— The observed oscillation periods can be formally divided into two groups: those 
with shorter (<2-min) and those with longer (>2-min) periods. 

— The most frequently observed oscillations are with period ranges of 50 — 110 s and 
3 — 7 min. 

— The propagation of the actual oscillations is rather difficult to detect. However, the 
relative Fourier phase between oscillations at different heights indicates the prop- 
agation speed of ~ 110 km . In some oscillations, perhaps caused by standing 
patterns, waves seem to be present with very high phase speeds (> 300 km s~^). 

The observed oscillations in spicules are most likely due to the propagation of 
transverse waves from the photosphere towards the corona. There are several possible 
interpretations of these oscillatory phenomena: 

— Kink waves propagating along slender magnetic flux tubes, where spicules are 
formed on fleld lines close to the axis. The kink waves lead to the transverse oscil- 
lations of spicules as a whole. 

— Volume-filling Alfven waves propagating in surroundings of spicules. The Alfven 
waves result in the oscillation of the ubiquitous magnetic field lines. These oscilla- 
tions force the spicules to periodic displacement of their axes. 

— Transverse pulses excited in the photospheric magnetic flux tube by means of buffet- 
ing of granules. The pulses may propagate upwards in the stratifled atmosphere and 
leave "wakes" behind, which oscillate at the kink cut-off frequency of the stratified 
vertical magnetic flux tube. The wakes can be responsible for the observed > 3-min 
period transversal oscillations in limb spicules. 

It is important to estimate the energy flux stored in spicule oscillations. If the 
oscillations are caused by the waves excited by granulation, then the energy flux of 
waves should be in the range of the energy flux of granulation. The energy flux stored 
in an initial transverse pulse at the photosphere is _F ~ UeC/^Vg, where Vg is the granular 
velocity, say 1-2 km s~^. Then, for the photospheric electron density of 2 ■10^'' cm^'^ 
and Alfven speed of 10 km s^^, the estimated energy flux is ~ 5T0^ erg cm^^ s^^ 
(taking Vg as 1.25 km/s). 

The same energy flux under spicule conditions (i.e. electron density of 2 -10^^ cm~'^ 
and Alfven speed of 100 km s~^) requires the wave velocity in spicules to be about 385 
km/s. The expansion of magnetic tube with height may alter the estimation. However, 
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spicule density is almost two magnitude higher compared to the other part of the tube. 
Therefore, even if spicules occupy only 1% of the tube cross section, the energy storied 
in spicule oscillations is comparable to the oscillation energy of remaining part of the 
tube. Therefore, even if the wave transmission efficiency is less than 100 %, the observed 
velocity is much lower than expected. This discrepancy can be resolved if the oscillations 
are caused by the wake behind a transversal pulse: in this case almost the entire energy 
of the initial perturbation is carried by the pulse, while the energy of the wake is much 
smaller. Therefore, even if the filling factor of magnetic tubes is 1%, the energy flux 
carried by pulses is more than enough to heat the solar chromosphere/corona. 



6.1 Targets for future observations 

More observations from space satellites and ground based coronagraphs are needed for 
a better and conclusive understanding of oscillatory events in solar limb spicules. There 
are few highlighted targets for future observations: 

6.1.1 Phase relations between oscillations in neighbouring spicules 

It is important to perform an analysis of phase relations between oscillations in neigh- 
bouring spicules. If spicule oscillations are caused by global Alfven waves or they are 
related to global photospheric oscillations, then the transverse displacement of spicules 
should show some spatial coherence i.e. a few neighbouring spicules should move in 
phase along the limb. Superposition of different spicule groups located along the line of 
sight could complicate the task, but careful analysis still may reveal some coherence. 
Hinode/SOT time series in Ca II H and Ha lines seem to be an excellent data for such 
analysis work. 

6.1.2 Phase relations between oscillations at different heights 

A study the phase relation of transverse displacements of a particular spicule at dif- 
ferent heights would allow us to infer whether the oscillations are due to standing or 
propagating wave patterns. Phase delays between different heights would determine 
the phase speed of perturbations, thus, the physical nature of the waves. Spectroscopic 
consecutive height series from ground based coronagraphs or time series of images from 
e.g. Hinode/SOT would allow to infer the wave length, phase speed or frequency of 
oscillations. These latter important diagnostic parame t ers m ay be then used to develop 
spicule seismology as suggested bv lZagarashvili et al. I l|2007l ). 

6.1.3 Propagation of transverse pulse 

Propagation of transverse pulses can be traced through a careful analysis of time series 
from e.g. Hinode/SOT. It is also important to search for oscillations of spicule axis at 
the kink or Alfven cut-off frequency as estimated in the chromosphere. Polar macro- 
spicules are probably the best targets for such work, as their life time is long enough 
when compared to the cut-off period. 
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